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Abstract 

The alpha-rich freezeout from equilibrium occurs during the core-collapse explosion of a massive 
star when the supernova shock wave passes through the Si-rich shell of the star. The nuclei are 
heated to high temperature and broken down into nucleons and a particles. These subsequently 
reassemble as the material expands and cools, thereby producing new heavy nuclei, including a 
number of important supernova observables. In this paper we introduce two web-based applications. 
The first displays the results of a reaction-rate sensitivity study of alpha-rich freezeout yields. The 
second allows the interested reader to run paramaterized explosive silicon burning calculations in 
which the user inputs his own parameters. These tools are intended to aid in the identification 
of nuclear reaction rates important for experimental study. We then analyze several iron-group 
isotopes (®®Ni, ^^Co, ^®Co, and ^^Fe) in terms of their roles as observables and examine the reaction 
rates that are important in their production. 
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I. INTRODUCTION 


Progress in the science of stars and their nucleosynthetic processes relies on the con¬ 
tinued interplay of astronomical observations and astrophysical modeling. Observations of 
abundances of chemical species, elemental or isotopic, constrain astrophysical models while 
the models, in turn, provide a framework for interpreting the observations. It has long 
been clear, however, that uncertainities in the input physics to the models limit them and 
their usefulness in interpreting abundance observations. Nuclear reaction rates are key in¬ 
puts into the astrophysical models, and, though many are measured, most have not, and 
modelers must therefore rely on theoretical predictions of the value of these rates. Recent 
theoretical reaction-rate predictions have proven fairly accurate (to within a factor of a few 
of the actual rate value where subsequently measured-e. g., |l|), nevertheless experimental 
results are usually desirable and often crucial. A third essential effort in nuclear astrophysics 
is that of nuclear experimentalists who seek to provide better input into astrophysical mod¬ 
els by measuring the rates of nuclear reactions or nuclear properties that improve theoretical 
estimates of the rates. Because of the cost in time, effort, and hnancial resources in perform¬ 
ing the necessary experiments, however, the task of the nuclear experimentalists is greatly 
aided if the astrophysical signihcance of a particular nuclear reaction can be clearly demon¬ 
strated. This in turn requires demonstration of the dependence of the predicted value of 
an astronomical observable in an astrophysical model on the value of the reaction rate. It 
is clear that observers, modelers, and experimentalists should all play an important role in 
the planning of nuclear astrophysics experiments, in particular by ensuring that any nuclear 
reaction rate proposed for experimental study satisfy the following four requirements P]: 

1. An appropriate astrophysical model of a nucleosynthesis process must exist. 

2. An observable from that process, usually an abundance result, is either known or 
measurable. 

3. The dependency of the value of the observable on the value of the nuclear cross section 
is demonstrable. 

4. An experimental strategy for measuring the reaction rate, or at least using experimen¬ 
tal data to better calculate the reaction rate, should be available. 
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A long-term goal of our program is to aid the dialogue among the three parties (observers, 
modelers, and experimentalists) in the nuclear astrophysics community by making it easier to 
identify those nuclear reactions most in need of experimental study. In the present work, we 
focus on reaction rates important in the alpha-rich freezeout in core-collapse supernovae. We 
explore in some detail astronomical observables from this process and the nuclear reactions 
that govern their nucleosynthesis. In addition, we present web-based tools that allow any 
interested researcher to explore reaction-rate sensitivities in the alpha-rich freezeout and 
thereby make a case for experiments on nuclear reactions important for other observables 
(yet to be identihed) from this process. We hope that this work can serve as a template for 
future work (by ourselves or others) on nuclear reactions important in other nucleosynthesis 
processes. 

This paper begins with an introduction to the alpha-rich freezeout and a description of our 
freezeout calculations and the attendant reaction-rate sensitivity studies in OT. In §pT| we 


introduce two web sites accessible from the main page jittp : //photon. phys. clemson. edu/| 
Ig jordan/nucleoTI . The hrst website displays data from the sensitivity calculation performed 
on our alpha-rich freezeout model in order to identify nuclear reactions that are important 
in the production of nuclei identihed as astrophysical observables. With a particular isotope 
in mind, one may then use the web site to view the effect of different reaction rates on the 
yield of that isotope. The second website can be used to calculate the effects of a varied 
reaction rate on nucleosynthesis yields under conditions that differ from the conditions used 
in our sensitivity survey. Control over several of the parameters in the explosive model is 
given to the user so that the effect of a particular reaction rate on the network yields can 
be explored over a wide variety of conditions, thereby allowing the user to strengthen his 
case for measuring a particular reaction. Finally, due to their signihcance as astrophysical 
observables, several isotopes from the iron-group nuclei are examined §|V[ 


II. THE CALCULATIONS 

Many processes contribute to the production of new nuclei in core-collapse supernovae, 
but one of the most important for astronomy is the alpha-rich freezeout from equilibrium. 
During the core collapse of a dying massive star, a shock wave develops as matter falls 
supersonically onto the collapsed stellar core. The shock, aided by a push from neutrinos 
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from the cooling nascent neutron star, expands out into, heats, and expels the overlying 
stellar matter. In the initial heating of the innermost regions of the ejecta, post-shock 
temperatures are sufficiently high that nuclei are broken down into nucleons and a particles. 
As the material subsequently expands and cools, the nucleons and a particles reassemble 
to form heavy nuclei. Because of the fast expansion of the matter, however, not all alpha 
particles reassemble, and, as a result, the hnal abundances freeze out with a signihcant 
number a particles remaining, hence the name alpha-rich freezeout. A number of signihcant 
astronomical observables are produced in this process including ^^Ti, ®®Co, and ^^Co. 

In order to explore the sensitivity of alpha-rich freezeout yields to variations in reaction 
rates, we utilized the Clemson nucleosynthesis code |Q], which we have updated to employ the 
NACRE 0] and NON-SMOKER 0 rate compilations. The network used (see [|]) includes 
376 species from neutrons up to Z = 35 (Bromine) and 2,125 reactions among them. 

We began the survey with calculations using reaction rates from the compilations. Guided 
by detailed models of such astrophysical settings (e. g., 0), we chose an initial tem¬ 

perature of Tg = T/10®K = 5.5 and initial density po = 10^ g/cm^ in all calculations. 
The matter was taken to expand exponentially so that the density evolved with time t as 
p{t) = po exp(—t/Tex), where the density e-folding timescale = 446 s/^po = 0.141 s. We 
assume the relation p (xT^ (radiation-dominated expansion), and took three possible values 
of the neutron excess rj, viz., 0, 0.002, and 0.006. Since we considered the appropriate model 
for alpha-rich freezeout to be the passage of a shock wave through Si-dominated matter, 
the inital composition was ^®Si with enough ^®Si to give the appropriate value of p. 

Upon completion of these reference calculations, we explored the signihcance of a partic¬ 
ular reaction by multiplying the reference value for its rate by a factor of ten, and repeating 
the calculation with the same expansion parameters as in the reference calculations. We 
emphasize that, with this procedure, we uniformly increase the rate at all temperatures 
by a factor of ten. The corresponding reverse reaction rate is also increased by the same 
factor at all temperatures. This is required by detailed balance to ensure the material 
reaches the appropriate equilibrium at high enough temperature and density. We repeated 
this procedure for all 2,125 reactions and all values of initial p’s. We also repeated this 
sequence with a multiplicative factor 0.1 on each of the reactions. The total number of 
alpha-rich freezeout calculations performed was (number of reactions) x (number of differ¬ 
ent values of p) x (number of multiplicative factors) -|- (number of reference calculations) 
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= 2,125 X 3 X 2 + 3 = 12,753. 

The effect of modifying a particular reaction can be determined from the ratio of the 
modified yield of a particular species with the reference yield. Clearly with 2,125 reactions 
and 376 species, there are many combinations to consider. Rather than present several large 
tables, we have opted to construct interactive web-based applications to display the results. 


III. INTERNET APPLICATIONS 

In this section we introduce two web-based applications. Since the web sites have complete 
online help files, detailed instructions are not given here. We will simply introduce the web 
sites and some of their relevant features, along with instructions on how to reproduce the 
tables appearing in the following sections. 

A. Network Sensitivity Data Display 


The Sensitivity Data Display website is designed to examine the results of varying the 
reaction rates involving a selected isotope. Upon loading the page, the user selects from 
one of six data sets. The data sets are specified by rj and by whether the reaction rates are 
increased or decreased (by a factor of 10 in each case). Next, the user can generate lists 
sorted by increasing or decreasing modi£ed-to-standard ratios for a selected isotope. An 
isotope (the chosen observable) is selected by entering its atomic number and mass number 
in the appropriate fields. A /3-decay option is also available which permits the user to evolve 
the abundances to arbitrarily later times upon cessation of the alpha-rich freezeout process. 
(The numerical technique used to construct the exact solutions to the coupled differential 
equations governing the /3-decay of the abundances is described in the Appendix). 

The data in Tables |V|, [V^, and |V11| of §|V| were collected using this web site. To recon¬ 
struct Table |VI1| , for example, the user selects the data set corresponding to r/=0.006 and a 
multiplicative factor of 0.1. For ®®Fe, 26 is entered as the Z value and 55 for the A value. 
We target an observation of the ^^Fe abundance made 2 years after the explosion, and so 
the user selects the /3-decay box and enters an evolution time of 2 years into the time field. 

In order to display only the reactions that change the ^^Fe yield by more than 20%, 
the ratio cut-off parameter is set to 0.20. The “Descending List” button is then clicked to 
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generate the ratio list for ^^Fe. The hrst table generated produced gives the mass fractions of 
all A = 55 species in the network at the end of the r] = 0.006 reference alpha-rich freezeout 
calculation together with their beta-decay or electron-capture lifetimes. Isotopes highlighted 
in black are stable while those in pink are unstable against or electron capture. Isotopes 
highlighted in blue are subject to /3“-decay while green denotes instability via both channels. 

In the second table, the column on the far right of the output contains the modihed- 
to-standard ratios after an interval of two years. The reactions which, when their rate is 
decreased by a factor of ten, result in a minimum increase in the yield of ^^Fe of 20% are 
entered in the table along with the respective ratios. Similarly, those reactions which, when 
their rate is decreased by a factor of ten, result in at least a 20% decrease in the ®^Fe yield 
were also entered in the table along with their ratios. The process was then repeated with 
a rate factor of 10 and rj of 0.006. The data for Tables and ^ were created in a similar 
manner, but for the isotopes ®^Co and ®®Ni. 

Other options on the web site include (a) a table of the isotopes used in network calcula¬ 
tions, (b) a list of reaction rates induced in the calculations and (c) a tool for the user to plot 
or make a table of the standard reaction rates used in the calculations. The user can use the 
reaction number from the reaction list and then input that reaction into the reaction-rate 
engine to obtain either a plot or table of the rate as a function of temperature. (Note that 
the rate we present may include a stellar enhancement factor to correct for the possibility of 
excited states in the target nucleus. The user must account for this factor when comparing 
his rate with the one we used.) Other links on the page are to various related topics in the 
help hie and to the NACRE and NON-SMOKER rate compilation web sites. 


B. Explosive Nuclear Burning 

The Explosive Nuclear Burning web site simulates parameterized explosive nuclear burn¬ 
ing over the web using the same reaction network as described above. Among the variable 
parameters are Tg (initial temperature in 10® K), p (initial density in g/cm^), r^x (the ex¬ 
pansion time scale of the explosion), and the initial mass fraction of each species. The user 
can also alter any reaction rate in the network by a multiplicative factor as was done in § ||. 

The data output from the web site are the hnal mass fractions from the calculation. If 
a reaction rate is modihed, the web site performs two calculations, one with the modihed 
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reaction rate, and one with the standard reaction rates for comparison. The data output 
in this case are the hnal mass fractions from both calculations along with their modihed- 
to-standard ratios. The web site offers several options for further analysis, including a (3- 
decay option to evolve the hnal mass fractions in the output. 

This web site was used to produce the data from Tables IQ.B and To recon¬ 
struct these tables, the user performs the following steps. First, the user selects the default 
parameters, which correspond to an alpha-rich freezeout of silicon-dominated matter with 
an Tj of 0.006 (the parameters the sensitivity calculations in § |T|). When the calculation is 
executed, the appropriate mass fractions are inserted into the “t = 0 years” column of the 
Table. The /3-decay box is then selected from the “Final Mass Fraction Data List Controls” 
table in the left hand frame and a time of 2 years entered. The “Create List” button is then 
clicked. The web site returns a table with the value of all of the hnal mass fractions after 2 
years. The hnal mass fractions should agree with those in Tables |, |1|, |T|, and |V. 

The data for Table |VIII was also generated with this web site. To reconstruct this table, 
the user performs the following steps. First, the web site is reloaded into the browser. Then 
the “Edit” button in the “Reaction Rates” row of the “Parameter Controls” table in the left 
hand frame is clicked. This brings up an interface that allows the user to select a reaction 
rate to modify based on the type of reaction preferred. Since the reaction of interest in 
Table |V llf| is ®®Co(p, 7 )®®Ni (a two nuclear species to one nuclear species reaction), the radio 
button at the bottom of the table entitled “z + j —^ k" is clicked. Starting at the left hand 
side of this table, the consituents of the reaction are chosen held by held. When all of the 
helds in the table are selected, the “Save Changes” button at top (or bottom) of the page 
is clicked. The user is then prompted to input a reaction rate multiplier. A value of 100 
is entered and the “Save Changes” button clicked. The web site then presents a summary 
of the changes and the calculation is begun by clicking the “Run Nucleosynthesis Code” 
button. The outputted data are then evolved for two years as above and the appropriate 


modihed-to-standard ratio is entered into Table [VIII| . The web site is then reloaded and 
the same procedure is carried out except 0.01 is entered as the multiplicative factor on the 
^®Co(p, 7 )®®Ni reaction rate. These steps should reproduce the data in Table IVB 

This web site has many other useful tools that are not described above but are documented 
in the on-line help hies. As shown, this web site allows the user to examine the effects of a 
reaction on nucleosynthesis yields under a wide variety of circumstances. Our hope is that 
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this easily accessible data will provide additional insight into the effects that uncertainties 
in reaction rates have on the yields of the alpha-rich freezeout and other nucleosynthesis 
processes. 

IV. NUCLEAR REACTIONS GOVERNING THE SYNTHESIS OF IRON-GROUP 
OBSERVABLES 

In this section, we present some results of our calculations drawn from the web sites 
described in §B For brevity we limit the discussion to the rj = 0.006 calculations. The 
purpose is to illustrate how one may identify relevant observables and identify their governing 
reactions using the web sites. 

At least three types of isotopic observables are relevant for stellar nucleosynthesis: 

1. The bulk yields are important for understanding galactic chemical evolution and solar 
system abundances. 

2. Radioactive species such as ^®A1 and ^^Ti can be observed from space telescopes, and 
consequently provide important constraints on their production sites. 

3. Isotopic abundances in presolar meteoritic grains carry important information about 
the nucleosynthesis environments in which their isotopes formed and the astrophysical 
settings in which the dust grains condensed (e. g.. I. i). 

In the subsections to follow, we present several isotopes from the iron-group nuclei for 
which we used the web-based applications in § 0 to obtain data. These four isotopes are 
directly or indirectly linked to all three types of observables listed above. We have chosen 
to analyze ®^Co and ^®Co because of their prominent role as y-ray observables and ^®Ni 
and ®®Fe for their potential as future X-ray observables. In addition, the abundances of the 
daughter isotopes of all four of these species may eventually be measured in presolar grains. 
We do not analyze ^“^Ti, a key alpha-rich freezeout observable, because this has been done 
previously p. 

For each species, we list any reaction rate that produces at least a 20% increase or decrease 
in the alpha-rich freezeout yield due to a factor of 10 change in the reaction rate (except for 
®^Co as explained below). Since we primarily consider the y-ray or X-ray observations of 


these species, we chose a time of analysis 2 years after the freezeout. By this time, overlying 
supernova should have become transparent to these radiations. Often, several years pass, if 
not thousands, before data from astronomical observables can be collected. It is, therefore, 
up to the user to make his own analyses using appropriately chosen times in the web tools 
described in §13 For each observable considered, we list the important reactions, and for 
two particularly interesting reactions we explore in some detail the reason for their effect on 
the chosen observable. 

A. 

This species is produced primarily as the radioactive parent ®®Ni (ri /2 = 6.075 days) 
which decays through ^®Co (ri /2 = 77.2 days) to ®®Fe (see Table |). The overlying supernova 
largely obscures y-rays from the decay of ^®Ni, but the 2.598 Mev y-rays (for example) from 
the ^®Co can escape and be detected [0, making the latter species a valuable alpha-rich 
freezeout observable||21||. The y-rays from ^®Co decay can be used to determine the total 
yield of the parent (^®Ni) as well as the opacities of the outer envelope of the supernova (for 
a review of y-ray observables see f^). 

From the web sites one hnds that the abundance of ^®Co after two years is dependent on 
the yield of ^®Ni from the supernova. We hnd that the yield of ^®Ni immediately after the rj 
= 0.006 alpha-rich freezeout is quite insensitive to factor-of-ten changes in the reaction rates. 
The largest effect was due to changes in the triple-a reaction which produced about a 6% 
change in the yield. This 6% change propagates as a 6% change in the yield of ®®Co after the 
decay of ®®Ni. These are not signihcant changes and we conclude that the calculated yields 
of ®®Ni (and thus the yields of ®®Co) are quite robust against any reaction rate uncertainties. 
This result is expected since the production of ^®Ni is dominated by the equilibrium phase 
of the expansion in which the abundances are set by binding energies and nuclear partition 
functions (not by individual reaction rates). 

B. 

A second important isotope is ®^Co. This isotope is similar to ®®Co in that its stable 
daughter ®^Fe owes a signihcant portion of its synthesis to alpha-rich freezeouts. From 
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Table H we see that is primarily produced as the parent isotope ^^Ni. Gamma-rays from 
the decay of ®^Co are observable with space detectors, and they also power the supernova 
light curve at somewhat later times than ^®Co |^. As shown, factors that affect the yield 
of ®^Ni thus have an effect on the yield of ®'^Co. 

The web sites show that the yield of ®^Ni immediately after the alpha-rich freezeout is 
(like ®®Ni) largely insensitive to the value of any particular reaction rate because both are 
primarily equilibrium products. The reaction with the most signihcant effect (for 77 = 0.006) 
is ^'^Ni(n,p)®'^Co. We have chosen to list this reaction even though it does not produce more 
than a 20% change because of the importance of ®^Co as an observable. From Table 0 we 
see that this reaction produces at most a 17.7% change in the yield of ^^Ni for a factor of 
ten change in the rate. The primary reaction rate governing the yield of ®^Co for r] = 0.006 
is ^^Ni(n,p)^^Co and a factor of a few uncertainty in the reaction rate results in a several 
percent uncertainty in the observable. 

The reason ®'^Ni(n,p)®'^Co has an effect on the yield of ^^Ni (and thus the yield of ®^Co) 
is that it governs when ®^Ni falls out of quasi-statistical equilibrium (QSE-see, for example, 
|r^). As shown in hgure |^, the ®^Ni abundance diverges from nuclear statistical equilibrium 
(NSE) early but remains in line with QSE expectations until it freezes out below Tg 
3.5. The QSE favors ^^Ni early, but as the material cools below Tg 4.3, the QSE abun¬ 
dances shift to higher-mass nuclei; thus, the longer ®^Ni remains in QSE, the lower its hnal 
abundance. ®^Ni falls out of QSE when the ^^Ni(n,p)^^Co reaction becomes too slow, which 
means that increasing the rate for this reaction causes ®^Ni to remain in QSE longer and to 
have a lower hnal abundance. Decreasing the rate means ^^Ni falls out of QSE earlier and 
retains more of its originally high QSE abundance. 


C. 59Ni 

Both ^®Ni and ^^Fe have been chosen because they hold promise as detectable radioactive 
X-ray sources (for a prospectus see [|^). In the decay of these isotopes by electron capture, 
a K-shell vacancy occurs which is often hlled by an electron previously occupying a higher- 
level bound state. The energy loss from this event is carried away as an X-ray which, if 
detected, would not only provide data complementary to 7 -ray observations, but also shed 
light on the production of radioactive isotopes that emit no 7 -rays. 
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The observable ®®Ni, because of its 75,000 y half-life, should produce a significant diffuse 
emission from the interstellar medium similar to that observed from interstellar ^®A1. It 
may also be detectable in close, individual supernova remnants |^. From Table H we see 
that the majority of ®®Ni is produced in the alpha-rich freezeout as ^®Cu; thus, factors that 
affect the production of ®®Cu ultimately affect the production of the observable ^®Ni. From 
the web sites it is evident that there are several reactions which, when changed by an order 
of magnitude, have a moderate effect on the production of ®®Cu (see Table [V^) . The most 
influential of these is the ®®Cu(p, 7 )®‘^Zn reaction. A factor of 10 decrease in this reaction 
rate produces about a 50% increase in the yield, while a factor of 10 increase in the reaction 
rate drops the yield of ®®Cu by about 30%. These are modestly significant changes and show 
that ®®Ni is sensitive to specihc reactions in the alpha-rich freezeout. 


D. 


55 


Fe 


Unlike ^®Ni, ^®Fe has a relatively short half life (ri /2 = 2.73 yrs), making it only detectable 
in very young supernovae (supernova 1987A remains a good target for detecting ®®Fe [Q). 
Since the detection of ®®Fe could become available at very early times after the explosion 
(assuming the ejecta is transparent to X-rays), its detection should give information on 
small-scale structure formation in the ejecta. Further, it may help elucidate the velocity 
structure of the core and the velocity structure of ®®Fe . 

As can be seen from Table ®®Fe is predominantly produced as ®®Co. The production 
of ®®Co in the alpha-rich freezeout with rj = 0.006 is most sensitive to the triple-a reaction 
and the ®®Co(p, 7 )^®Ni reaction (Table [VII|) . The triple-a reaction has about a factor of two 
effect on the yield when a factor of 10 change in the triple-a rate is made. The largest effect 
on the yield is seen in the dependence of the ®®Co abundance on the ®®Co(p, 7 )®®Ni reaction. 
A factor of 10 decrease in this reaction rate increases the yield of ®^Co by 526% while an 
increase in the reaction rate reduces the yield by almost 80%. 

The reason for the sensitivity of ^®Co to the ®®Co(p, 7 )®®Ni reaction rate is similar to that 
of the sensitivity of ®^Ni to the ®^Ni(n,p)®^Co reaction rate, but the effect is much larger. 
The ^®Co remains in (p,7)-(7,p) equilibrium with abundant ®®Ni until Tg drops below 3.0. 
As the temperature drops, the equilibrium shifts abundance from ®®Co to higher-mass nuclei 
(particularly ®®Ni). Increasing the ®®Co(p, 7 )®®Ni reaction rate causes this abundance shift to 
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persist longer, which in turn leads to a lower final ^^Co yield. A lower value for the reaction 
rate leads to an earlier freezeout from QSE and a higher final ^®Co yield. From this we see 
that ®®Co, and thus ^^Fe (the relevant observable), are very sensitive to the ^^Co(p, 7 )®®Ni 
reaction rate. 


V. CONCLUSION 


The above described web sites are intended to be a starting point for identifying those 
nuclear reactions that govern the production of astrophysical observables from the alpha-rich 
freezeout. It is our hope that these web sites will help motivate future measurements of these 
key reactions, especially as new observables from the alpha-rich freezeout are identified. 

While we identified and analyzed four key observables from the alpha-rich freezeout, 
we expect new ideas or emphases in astrophysics to make new observables available. For 
example, as models of Galactic chemical evolution become increasingly refined, better yields 
from massive stars will become necessary. This means that any isotope that owes a significant 
portion of its solar system abundance to the alpha-rich freezeout, such as ^“^Ca, ^®Ti, 52,53 Qj, 
and the ones discussed above, will become an observable. Similarly, new technologies may 
also open new possibilities for observing isotopes. As a possible example, ^^Ca excesses 
found in presolar supernova silicon carbide (SiC) X grains arose from condensation of live 
^^Ti [^. One speculation is that much of this excess radiogenic ^^Ca may be concentrated in 
small titanium carbide (TiC) subgrains within the larger SiC X grains. Such TiC subgrains 
are known to exist in mainstream SiC grains that condensed in outflows from low-mass 
stars 1^, however, in two SiC X grains studied so far, there is no direct evidence for 
the presence of such TiC subgrains[^], though the search continues. If such subgrains 
do exist, they are likely to be dominantly comprised of alpha-rich freezeout material. If 
the new generation of secondary-ion mass spectrometers are able to measure the isotopic 
abundances in the subgrains, then the alpha-rich freezeout abundances of all titanium and 
carbon isotopes, as well as the isotopes of any element that might condense in TiC, would 
become observables. Those new observables will open up new governing reactions requiring 
experimental study. 

We view the web sites presented here as a first step for similar, future work on other 
nucleosynthesis processes. The goal in those future efforts, as in the present work, will be 
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to help identify nuclear reactions that govern the nucleosynthesis of astronomical observ¬ 
ables. As always, however, the challenge remains not only to identify nuclear reactions that 
govern the synthesis of those observables, but also to identify the key astronomical observ¬ 
ables themselves. Only by maintaining a healthy dialogue among astronomers, astrophysical 
modelers, and nuclear experimentalists can we face that challenge and advance the science 
of nuclear astrophysics. 
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APPENDIX: ANALYSIS OF THE /3-DECAY ALGORITHM 

We now draw attention to the algorithm used by the Sensitivity Data Display and the 
Explosive Nucleosynthesis web-sites to compute alpha-rich freezeout yields at a specihed 
point in time. Traditional approaches rely on Runge-Kutta or implicit Euler techniques to 
integrate the coupled differential equations (preferably with adaptive timestep adjustments). 
However, we have sought a rapid, but exact solution. To do this, we exploit the elementary 
matrix solution to the system of linear coupled DEs that utilizes eigenvalues and eigenvectors. 
For any matrix equation of the form 


X = Ax 


(A.l) 

(A.2) 


x(0) = b 


the solution at all later times is given by: 



(A.3) 


14 


where are the eigenvectors and ook the corresponding eigenvalues of the matrix A 
(constructed from /?—decay rates and the application of mass conservation) which remain 
constant during the evolution. The coefficients Ck of the linear combination are the elements 
of the product (where T is the matrix of eigenvectors of A ordered strictly according 

to the index “fc”) 

The technical challenge here is that A is singular and thus reduction to upper Hessenberg 
form (the practical method of computing the spectrum of a large asymmetric system) using 
QR Factorization fails. We may circumvent this difficulty by employing very small perturba¬ 
tions to the rates and thus disturbing the rate matrix away from singularity, as is the norm 
for practical Numerical Linear Algebra applications in most fields of engineering. However, 
we must be cognizant of the fact that a theorem of Watkins requires the eigenvalues 
to be reasonably separated for the system to be insensitive to perturbations. Unfortunately 
for the full 376x376 matrix, the observed degeneracies are quite severe. 

On closer inspection, however, it is realized that the cause of the degeneracies is quite 
simple: the matrix is composed of several non-communicating blocks. Since only /?—decays 
are considered, the vector space decomposes into the direct sum of these blocks, each cor¬ 
responding to a particular mass number. Blocks comprised of isotopes from different mass 
numbers cannot traffic with each other since we exclude other kinds of nuclear reactions. 
Several of these component submatrices are very similar in their spectra, and therefore the 
full space suffers from many closely spaced eigenvalues. 

The solution, however, is now obvious-we merely compute the spectrum of each block, 
confident in the knowledge that the isotopes within a particular block are not affected by 
those in any other. The complete state vector of the abundances at any given time is then 
a simple concatenation of the state vectors from the subspaces. Watkins’ theorem is now 
observed to hold for each block (though the blocks themselves are singular as the full matrix 
was before) and abundances are output within seconds. 

Obviously, the next step is to extend the perturbation technique to matrices which can¬ 
not be decomposed into non-communicating subspaces (blocks). Weakly coupled blocks may 
lend themselves to similar perturbation analyses. Unfortunately Watkins’ theorem does not 
predict error bounds and the most we can do at this stage is carry out some rough numerical 
experiments to check the sensitivity of the system to small random perturbations. In our 
trials we found the effects to be negligible. Propagation of the perturbations caused fluctu- 


15 



ations on the order of 10“^^ in the most ill-conditioned blocks. Abundances of order 10“^® 
or less are quite small and for any practical purpose such small abundances are nonexistant. 
We conclude that the system was very robust to the small kicks we gave it. 

We hope to incorporate the latest techniques in spectral computations on large asymmet¬ 
ric systems as they become available, since nucleosynthesis would be an exciting test-bed 
for these tools from the frontlines of research in Numerical Linear Algebra. In the interim, 
perturbation allows us to practically carry out an important calculation in seconds which 
would otherwise have been impossible if we were to insist on a spectral solution. 
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TABLE I: Mass Fractions for A=56 with rj = 0.006 


Isotope 

t = 0 years 


t 

= 2 years 

56Ni 

0.767 



0 

56Co 

3.356 X 10-5 



1.190 X 10-5 

56Fe 

5.161 X 10-12 



0.766 


TABLE II: Mass Fractions for A=57 with r/ 

= 0.006 



Isotope 

t = 0 years 


t 

= 2 years 

57CU 

0 



0 

57Ni 

3.967 X 10-2 



0 

57Co 

6.735 X lO-"^ 



6.197 X 10-5 

57Fe 

3.429 X 10-1^ 



3.348 X 10-2 


TABLE III: Mass Fractions for A=59 with r] 

= 0.006 



Isotope 

t = 0 years 


t 

= 2 years 

59Ga 

0 



0 

Zn 

6.260 X 10-1^ 



0 

59Cu 

1.887 X 10-5 



0 

59Ni 

4.065 X 10-5 



1.928 X 10-5 

59Co 

3.682 X 10-15 



3.562 X 10-5 


TABLE IV: Mass Fractions for A=55 with rj 

= 0.006 



Isotope 

t = 0 years 


t 

= 2 years 

55Ni 

3.306 X 10-15 



0 

55 Co 

1.013 X 10-5 



0 

55Fe 

3.393 X 10-15 



6.058 X 10-5 

55Mn 

0 



4.072 X 10-5 
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TABLE V: Data from the rj = 0.006 network survey for ®'^Co 


Reaction 

Reaction Rate x 0.1 

Reaction Rate x 10.0 

®'^Ni(n,p)^'^Co 

1.177 

0.876 


TABLE VI: Data from the rj = 0.006 network survey for ^®Ni 

Reaction 

Reaction Rate x 0.1 

Reaction Rate x 10.0 

59Cu(p,7)‘^0Zn 

1.556 


0.615 

2a{a)^^C 

1.324 


0.731 

^®Cu(p,Q;)^®Ni 

0.848 


1.444 


TABLE VII: Data from the rj = 0.006 network survey for ^^Fe 


Reaction 

Reaction Rate x 0.1 

Reaction Rate x 10.0 

55Co(p, 7) ^‘^Ni 

5.258 

0.232 

59Cu(p,a)56Ni 

1.170 

0.738 

^®Cu(p,7)®'^Zn 

0.690 

1.133 

2a(Q;)^^C 

0.561 

1.887 


TABLE VIII: Data from web based Explosive Silicon Burning with rj = 0.006 for ^^Fe 


Reaction 

Reaction Rate x 0.01 

Reaction Rate x 100.0 

55Co(p,7)56Ni 

34.544 

5.469 X 10-2 
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FIG. 1; Mass fraction of ^'^Ni versus Tg for different values of the reaction rate ^'^Ni(n,p)®'^Co and 
in NSE and QSE. The value of the reaction rate near Tg ~ 3.5 governs when ^^Ni breaks out of 








FIG. 2: Mass fraction of ^®Co versus Tg for different values of the reaction rate ^^Co(p, 7 )®®Ni and 
in (p,7)-(7,p) equilibrium. The value of the reaction rate near Tg ~ 2.5 determines when ®^Co 
breaks out of (p,7)-(7,p) equilibrium with ^®Ni. 







